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Abstract 

We point out that the disagreement between the paper by Neukirch and Starostin (Ref.[l]) and ours 
(Ref. [5]) is only apparent and stems from a difference in approach. Ref. [1] is concerned with classical 
elasticity and individual curves while Ref. [5] focuses on statistical averages over curves. 
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Recently Neukirch and Starostin [l| have noted that the use of Fuller's formula for writhe may 
not always be justified in analysing experiments that stretch and twist single DNA molecules. 
They have criticised some earlier works [17-34 of Ref. for using Fuller's formula sl, 4] without 
always being careful to check the conditions for its validity. Readers of this paper [l] may derive the 
impression from some remarks made about 5], that Ref.[l| somehow invalidates the conclusions 
of our paper j^, which shows how Fuller's formula can be used in understanding entropic DNA 
elasticity. Here we note that the apparent differences can be traced to a difference in approach 
to the problem. We will set the comment in perspective by explaining the two points of view 
that have influenced the literature on bio-polymer elasticity. We summarise the discussion by 
concluding that the remarks made in [l| in no way invalidate any of the claims made in j^. 

Two communities with slightly different approaches have been working on bio-polymer elastic- 
ity. One point of view (See, for instance, Q]) is purely mechanistic and has its roots in classical 
elasticity. While this point of view captures some of the qualitative features of single molecule 
experiments and works well in the energy dominated regime of stiffer polymers like Actin and Mi- 
crotubules, it fails to capture the regime where there is a competition between the intrinsic elastic 
energy of the polymer and thermal fluctuations which are present in a real cellular environment. 
This latter regime which is better explained from a statistical mechanical point of view (See, for 

Q 

instance, [7|), where the central notion is the partition function of the system. These two com- 
munities view bio-polymers from slightly different perspectives: the mechanistic view emphasizes 
individual configurations, while the statistical view averages over configurations and focuses on 
the partition function, which is related to experimentally accessible quantities. The differences 
between our paper (Ref [5]) and Neukirch and Starostin's paper (Ref[l]) comes from these two 
distinct viewpoints. 

The view offered in (Sfl, is that while it is incorrect to claim that the partition function of the 

n 

SAWLC model equals that of the worm like rod chain WLRC exactly (as [7[ appear to do), it 
is nevertheless a good approximation over a range of forces and torques. To see this, note that 
changes in the two notions of writhe (Fuller Writhe Wp and Calugareanu- White writhe Wqw) are 
equal {5WcwiC) = 6Wf{C)) for small variations of the curve C (provided both quantities are well 
defined in the variation). Integrating this equation we find that the difference Wcw{C) — Wf{C) 
is constant for deformations of C which are neither self crossing nor south crossing. We will follow 
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[l| in referring to such deformations as "good" deformations. We choose for a reference curve the 
straight hne in the z direction it = z). Noting that the constant vanishes on the reference curve, we 
arrive at Fuller's formula Wp = Wcw for all curves which can be deformed to the reference curve 
by "good" deformations. We refer to curves related to the reference curve by "good" deformations 
as "good" curves. Note that the set of "good" curves is much larger than just small perturbations 
about the straight line. For instance, curves which are nowhere back-bending {tz > 0) are "good" 
curves and these may be far from straight. 

The main point we make in Ref [5] is that for a range of {F,W), the set of "good" curves 
dominates both partition functions and as a result, Zsawlc{Fi W) equals Zwlrc{F, W) approxi- 
mately. 

ZsAWLciF, W) ^ ZwLRciF, W) (1) 

The accuracy of the approximation is determined by the extent to which the "good" curves dom- 
inate the partition function. A simple example clarifies the matter. Suppose we wish to find the 
expectation of the function h{x): h{x) = x"^ when —3 < a; < 3; h{x) = — 1 otherwise, defined 
piecewise over the real line with a Gaussian measure exp{—ax^), suitably normalised. Consider 
the function g{x) = x"^ everywhere. Certainly h{x) is not equal to g{x), as there are points where 
they differ considerably. However, in computing the expectation value < h{x) > of h{x) one 
can approximately replace it by < g{x) > which is analytically more tractable. < h{x) > is 
approximately equal to < g{x) >. The approximation is good if a is not too small. Likewise, for 
forces which are not too small, the approximate partition function is expected to be close to the 
exact one. This explains the efficacy of Bouchiat and Mezard's WLRC model beyond perturbation 
theory y,y. 

n n 

Both Ref.pj and Ref.[5| attempt to understand the two writhe formulae and their use in the 
non perturbative regime. 

n 

However, the emphasis in p| is on individual configurations and not in the statistical sense. 
It is indeed true as both Ref. [1] and Ref. [5] note, that the writhe formulae are not the same 
on individual curves. Our claims of approximate equality of the two models are made at the 
level of the partition function and not individual curves. We claim that the partition functions 
are approximately equal (Eq. ([1])) to each other for a wider range of parameters than one would 
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naively expect. Such an approximation works very well at very high forces. Let us now consider 
twisting the molecule: 

(a) . At low twist: This is the paraxial limit where the backbone of the polymer is essentially 
straight and the tangent vector i to the polymer explores the neighborhood of the north pole 8|, l9| . 
This is the perturbative regime, which is not controversial. 

(b) . At intermediate twist: In this regime we have a writhing polymer which may not be nearly 
straight. However, our argument that "good" curves dominate the partition function applies and 
we conclude that Eq. ([T]) holds. 

(c) . At high twist: The energy cost of accommodating writhe becomes nearly zero in self 
avoiding models. This is because an infinitely thin polymer can writhe at negligible energy cost 
by winding around itself as a plectoneme. (The word is Greek for "twisted thread" and describes 
the structures often seen on telephone cords.) To render the energy finite, one has to "fatten" 
the thread and allow for the finite thickness of the DNA molecule (about 2nm). This pathology 
of infinitely thin threads is well known to mathematicians and was quite early noticed by Fuller 

A similar pathology also afflicts south avoiding models: writhe can be stored at negligible 
energy cost by winding around the south pole. These configurations can be described, mixing our 
small Latin and less Greek, as "Australonemes" (southerly threads). If one excludes a finite region 
around the south pole (as Bouchiat and Mezard do 3]) by using a cutoff, one ends up with the 
same finite energy cost per unit writhe for appropriate cutoff. As a result, even in the high twist 
limit one finds that the partition functions are approximately equal and Eq. ([1]) holds. 

To summarize, we have clarified the issues surrounding the use of Fuller's formula for the 
writhe of a space curve. We note that the formula, when used with due care, can be a valuable aid 
to taming an otherwise intractable calculation. It permits an approximate determination of the 
partition function and a prediction for the experimentally measured twist-extension relations in the 
presence of an applied force. We observe 5|] using the Closed Circuit theorem that self avoidance 
and south avoidance have the same topological effect of obstructing link release and topological 
untwisting. Our observations justify theoretical work by Bouchiat and Mezard, which though 
successful in interpreting the data, has been criticised for the incorrect use of Fuller's formula. 
We thus understand the "unreasonable effectiveness" of Fuller's formula in understanding DNA 



4 



elasticity. 
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